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Abstract. We study the irreducible and algebraic equation 

(1) x"+4a 1 x"- 1 + ha„=0,n>4 

on a differential field (F = C(x),8 =') of characteristic zero with algebraically closed field of 
constants C. We assume that a root of the latter equation is solution of a Riccati differential 
equation 

u' = B + Biu + B 2 u 2 

where Bo, Bi, B2 are in F. 

We obtain a parameterization of the polynomials in terms of a variable T g F. All the possible 
solutions are of the form 

F(x,T) = 

with numerical coefficients. We also give the possible algebraic Galois groups in each case and 
build explicitly the concerned Riccati equations. We then give an example of a degree 3 ir- 
reducible polynomial equation satisfied by certain weight 2 modular forms (for the subgroup 
T(2) = {M g SL(2, Z), M = Id mod 2} of 51/(2, Z)), all whose solutions satisfy a same Riccati 

equation on the differential field ( C(£?2, E4, Eg), — ) with Ei the Eisenstein series of weight i 

V drj 

respectively. These latter solutions are related to a Darboux-Halphen system. We then show 
that in the generic situation, corresponding to when the three solutions of the Darboux-Halphen 
system are all distinct, these latter satisfy an irreducible degree three polynomial on the dif- 
ferential field \C(-y, 7', 7"), — J , with 7 the generic solution of the Chazy differential equation 
\ dry 

7"' = 67"7 - 9-y' 2 . Finally our interest is the following problem: for which potential q does the 

Riccati equation \- v? = q admit algebraic solutions over the differential field C(p, p'). In 

dz 

this problem, the approach will be through Darboux polynomials. 
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1. Introduction 

Let (F = C(x),5 —') be a differential field of characteristic zero with algebraically closed field 
of constants C. We proceed in this paper to the study of some particular properties of the Riccati 
equation which is a differential equation of the form signaled above. We consider an algebraic 
irreducible equation h n (x) of degree n > 4 ([1} 

(2) h n {x) = x n + a^x^ 1 + ■ ■ ■ + a n (t) = 0, n>4 

where the dj reside in the field F. Suppose a root of h n is solution of a Riccati equation 

(3) u'(t) = B Q + B lU + B 2 u 2 , 

where Bq, B\, B 2 are in F. As we are in characteristic zero, every irreducible polynomial is 
separable and hence has got a transitive Galois group H. So if one root of equation @ satisfies the 
equation all its roots do. Therefore our hypothesis is identical to giving oneself an irreducible 
separable h n (x) all whose roots satisfy a same Riccati equation of the form (j3|). 

The cross-ratio of any four solutions xi ([23]) will be a constant and we will call such equations, 
anharmonic equations. This fundamental property leads to profound implications; it is indeed 
this Constance which will enable us to effectively parameterize the coefficients of the anharmonic 
polynomials. In the following we will not distinguish between the h n © which differ only from a 
substitution 

fA . Ax + B 

(4) *^CZ+D 

with AC - BD ^ and A, B, C, D e F. 

This identification makes sense because such a transformation changes a solution of a Riccati 
equation to another solution of a Riccati differential equation and preserves the cross-ratio of any 
four elements in some field extension of F and does not change the irreducible nature of h n . 

We will take advantage of similar substitutions in order to assume that the sum a± of the roots 
of / is zero or to multiply x by an appropriate element of F. The goal being the simplification of 
the expressions by associating to a given anharmonic equation an equivalent simpler one. 

This is the main result that we obtain 

Theorem 1. Every h n is of the form 

F(x,T) = 0, 

where the polynomial in two variables F has constant coefficients (elements of C) and T 6 F. Its 
algebraic Galois group is a finite subgroup of PGL(2,C). 

The interest in the Riccati equation has always been huge; it is the only first order differential 
equation (solved in the first derivative) up to fractional linear transformations of the dependent 
variable and holomorphic transformation of the independent variable, possessing the Painleve prop- 
erty. To be more precise, one knows from the Cauchy theorem of existence and unicity of a local 
solution of a differential equation, that the general solution of a first order differential equation is 
of the form 

f(z,Co), with Co an arbitrary constant of C. 
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But in general the solution y(z, Co) exhibits singular points, meaning points where y(z, Co) is not 
analytical. These singularities are of different types: poles, or possibly less friendly ones called 
branch points and essential singularities. 

In most cases the location of these singularities depends on the constant Co- The Riccati 
equation is the only first order equation whose only movable singularities are poles (Painleve 
property). One may consult [1J| for the latter subject. 

Another aspect of the theory of Riccati equations is the linearization property. In fact the 
substitution 

z' 

changes the Riccati equation into a second order linear one. Therefore the particular Riccati 
equation 

is related in this way to the one dimensional Schrodinger equation. 

d 2 z 

Nowadays the Riccati equation is most notably used in the field of linear differential Galois 
theory where it is used in order to characterize existence of Liouvillian solutions for the second 
order linear differential equation This paper is structured as follows: in section[5]we give 

a review of differential Galois theory, prove that the Galois groups of anharmonics are subgroups 
of PGL(2,C), give the explicit form of the Riccati equations and end with the main theorem. Then 
we give the possible degrees of anharmonics and example in subsection 12.11 Then we study the 
irreducible degree 3 polynomials coming from modular forms and the Darboux-Halphcn system. 
Finally in the last section [5] we are interested in the following problem: for which potentials q does 
dy 

the Riccati equation - — h y 2 = q admit algebraic solutions over the differential field C(p, p'). 
dx 



2. Algebraic solutions of the Riccati equation 

We first remind some facts about the Galois theory of second order linear differential equations. 
One may consult [27| for a more complete exposition about the subject. 

A differential field (F, 5) is a field together with a derivation <5 on F. One says that a differential 
field (J, A) is a differential field extension of (F, 5) if J is a field extension of F and its derivation 
A extends S. 

In the following y' and y" will stand for 8(y) and S 2 (y). Let an ordinary homogeneous second 
order linear differential equation 

L{y) = y" + b lV ' + b y = 0, o,eF. 

Let ?y, C a fundamental set of solutions of L (two independent solutions generating its two dimen- 
sional vector space of solutions V over C). Form the differential extension field 

3 = W <7 7 ,C>=F(r 7 ,? ? / ,C,0- 

The former field extension $ is called a Picard-Vessiot extension if in addition 5 and F have the 
same field of constants. The differential Galois group <5(L) of 5 over F is the group of differential 
automorphisms of $ that leave F invariant. (An automorphism a is differential if a (a 1 ) = {era) 1 , 
a €3)- 



4 



AHMED SEBBAR AND OUMAR WONE 



The previous choice of a basis of solution induces a faithful representation of <8(L) as a subgroup 
of GL(2,C), defined in the following way: for u E &{L) one has 

\a(ri) = a a r\ + b a ( 
\<t(Q = c a 7] + d a (. 

A different different choice of basis (171, Ci) leads to equivalent representations as there exists 
M € GL(2,C): (771, Ci) — M(r),Q. We identify these equivalent representations. Also one can 
show, (see [271 . p. 19]), that the differential Galois group considered as a subgroup of GL(2,C) is an 
algebraic subgroup. 

We have the following Galois correspondence [27, p. 25] 

Proposition 2. Let L(y) = y" + a\y' + a 2 y — be a second order equation over F with Picard- 
Vessiot field $ and &(L) :— Gal($/¥) its differential Galois group. Consider the two sets 

& := the closed algebraic subgroups o/©(L) (in the Zariski Topology) and 
£ := the differential subfields M of containing F. 
Define a : & — > $ by a{H) = $ H , the subfield of H -invariant elements of$. Let f3 :£,—>■ & defined 
by (3{M) = Gal($/M), be the group of M -linear differential automorphisms. Then 

(1) The maps a and f3 are mutual inverses. 

(2) The subgroup H £ & is normal if and only if M — $ H is, as a set, invariant under G. If 
H € & is normal then the canonical map ©(L) — > Gal(M /¥) is surjective and has kernel 
H . Moreover, M is a Picard-Vessiot field for some linear differential equation over F. 

(3) Let <5(L)° denote the identity component of <&(L). Then D ¥ is a finite Galois 
extension with Galois group <8(L)/<5{L)° and is the algeb raic closure o/F in 

y' 

Now the mapping y M- — defines a surjection between the set of non trivial solutions of L(y) = 

y 

y" + aiy' + a2y — and the solutions of the Riccati equation u' + b\u 2 + b^u + 63 = 0, with values bi 
function of and their derivatives. In the sequel we assume we deal with this form of the Riccati 
equation, which we suppose to have solutions which are algebraic of degree> 4. We remind that 
h n {x) = x n + aix n ~ 1 + • • • + a„, e F denotes its minimal polynomial. Let Xi, i = 1, • • • , n be its 

y' y' 

roots in a decomposition field. There exists — , i = 1, • • • , n: Xi = — . Using a fundamental basis 

y% yi 

of solution (77, C) of y" + a\y' + a-iy, one sees that the Ui all belong to the Picard-Vessiot extension 
"$ = ¥ < rj, C >= ¥(r),ri' ,(,(,'). But the field extension ¥[xi, i = 1, • • • , n] of F is a differential 
subfield of the Picard-Vessiot extension as x' i + b^x\ + b 2 Xi + 63 = 0. Thus it corresponds by 
the Galois correspondence to a finite closed algebraic subgroup H of the differential Galois group 
&(L) C GL(2,C). These differential automorphism of ¥[xi, i = ,n] are in the algebraic 

Galois group of h n . Conversely let a belong to the algebraic Galois group of h n ; by hypothesis it 
is transitive and permutes the roots xt of h n . If cr(xi) = Xj for i 7^ j, the case i = j being clear, 
one has 

a(xi)' = x'j = - (hx? + b 2 Xj + b 3 ) 

and 

a(x'i) = a(-(bixf + b 2 Xi + 63)) = -bia(xi) 2 - 6 2 cr(a; i ) - b 3 = -b\x 2 - b 2 x 3 - 63; 

therefore extending to the algebra F[a;,, i = 1, • • • ,n], one sees that any element of the algebraic 
Galois group is also differential and fixes F, ie belongs to H . Now one remarks that any element 
of &(L) the subgroup of those a € ®(L) that act on the solution space V as scalar multipli- 

y' 

cation, has a trivial action on any — , for y any non-trivial solution of L(y) = y" + a\y' + a 2 y = 0. 
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Thus one has, passing to quotient space the following lemma (as C is algebraically closed one 
has PSL{2,C) = PGL{2,C)). 

Proposition 3. The Galois group ofanyh n is a finite algebraic subgroup of PSL{2,C) = PGL{2,C), 
with algebraic meaning that it is the image under the quotient map of a finite algebraic subgroup 
ofSL(2,C). 



The list of the finite algebraic subgroups of SL(2,C) is furnished in [171. Il8|. We have 

Theorem 4 ( E3| ) • Let G be a finite subgroup of SX(2,C). Then either 

• G is conjugate to a subgroup of the group 

with D the diagonal group, or 

• the order of G is 24 {the tetrahedral case), or 

• the order of G is 48 {the octahedral case), or 

• the order of G is 120 {the icosahedral case). 

In the last three cases G contains the scalar matrix —1. 

These groups admit well-known representations as fractional linear transformations (see fl7|). 
The tetrahedral group A4 of PSL{2, C) is isomorphic to the group of linear fractional transforma- 
tions generated (under composition) by the elements 

B 2 : z -> -z; e : z -> -; &i : z, — — - 

z 1+jz-I 

the dihedral group 2) m is isomorphic with the subgroup of linear fractional transformations gen- 
erated by 

9 m : z — > £z; eo : z — > -, where ^ m = 1. 

z 

The octahedral group is on his side generated by the two linear fractional transformations 

1-iz+l 1-t 

Sl'.Z-^——. t, 2 :z^——z; 

1 + i z — 1 1 + 1 

the icosahedral comes from the linear fractional transformations 

B5 : z — > £,z Co : z — > — , e\ : z — > — , with f = 1, a = =— and p — -= — . 

z pz-a V5 V5 

Finally the cyclic subgroups of order n are generated by 

0„ : z -J- C = 1. 

To each of the previous groups of fractional linear transformations is associated an absolute invari- 
ant or "automorphic form". 

For the tetrahedral case, an absolute invariant is given by the formula 

, ,_rl>(z) _ f z± + l + 2iV3z 2 \ 3 
^ Z) -W)-\z^l-2^j ■ 

For the dihedral cases one can consider the map 

j \ ip {z) 1 _ 
4>{z) z m 
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Absolute invariants for the octahedral and icosahedral groups read respectively 

V; (g) _ (1 + 14z 4 + z 8 ) 3 
oct[Z) 4>(z) I08z 4 (I-z 4 ) 4 

and 

jj(z) _ (-(z 20 + 1) + 228(z 15 - z 5 ) - 494z 10 ) 3 
lco[z) ~ <f>(z)~ I728z 5 (z 10 + Hz 5 - l) 5 

The cyclic group admits the invariant 

Vc?/clzj " d>(z) 

The properties of the absolute invariants are the following: 

• They are left invariant under the transformation of the concerned groups, with the action 
given by 

• Every other absolute invariant for the group is a rational function with coefficients in C of 

Before continuing we state the Constance of the cross-ratio of four solutions of the Riccati equation 

Lemma 5. Let (F, 5 =') a differential field and u' + b\u 2 + b 2 u + b 3 = a Riccati differential 
equation with coefficients in F. Let U\, u 2 , u 3 , u 4 four arbitrary distinct solutions of the Riccati 
equation. Form the cross-ratio 

(ui - u 3 ){ui - u 4 ) 
(U2 - u 3 )(u 2 - Ui) 
and take the derivative with respect to 8; then it vanishes identically. 

These absolute invariants will be useful in order to express the roots of the anharmonics. We 
will then use the Constance of the cross-ratio to determine a form for the general solution of the 
Riccati equation and finally build the Riccati equation into consideration explicitly. 

For each anharmonic h n with algebraic Galois group H C PSL(2,C) (as represented above), 
consider one of its roots, x i7 which is therefore solution of the associated Riccati equation u' + 
b\u 2 + b 2 u + 6 3 = 0. Let H C PSL(2, C) the Galois group of h n and t a variable then 

*(t) 

is left invariant by H; therefore belongs to F. Consider Te F \ C. 

Let H Xi the stabilizer of a root of h n . As the Galois group of h n is transitive, one has if N 
denotes its order, N = np, where p is the cardinal of the stabilizer. Consider the action of H Xi on 
the set of roots Xi of h n . Two cases can happen. H Xi does not fix any other root xj 7^ Xf, then 
its action on a root Xj ^ Xi gives rise to an orbit H Xi Xj of length p. By the disjointness of the 
different orbits and the fact that the total sum of their cardinals is n — 1, one has: p\(n — 1). In 
the same way one sees that p\(n — 2) if it happens that H Xi fixes another root Xj ^ Xi. 

No element of H different from the identity element can fix more than two roots by Constance 
of the cross-ratio. Regarding now H as an homography group, we denote by i = 1, • • • , n a 
fixed-point of the corresponding homography subgroup of H Xi , Si. We first have the following 
lemma which shows existence of fixed points for finite homography subgroups. 



Lemma 6. Every finite homography subgroup admits a fixed point. 
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Proof. Let the equation 



= x, ad — be 7^ 



hold; then one has 



cx + d 
cx 2 + (d — a)x — b = 0. 



If c 7^ then a fixed point exists as the field of constant is algebraically closed. If c = 0, then if 
d =/= a a fixed point exists. If d = a ^ 0, then the group is infinite unless 6 = 0, which correspond 
to the identity (the field is of characteristic zero). Hence the result. □ 



The irreducible equation 

(5) 

has got 



*(t) 



m 
0(t) 



T, 



N 



np 



roots in an algebraic closure. We suppose that their list is given by t,-, where < j < N 
Consider the equation 



1. 



(6) 



i>(v) V'Oto) ip(vi) 



Hv) <t>(vo) 4>{m) Hvn-i)' 

If rjj is one root it, then for any a in its stabilizer Sj, r\j is a solution of the equation 



(7) 



<t>(a.r)) 4}(rjj)' 



therefore rjj is root of equation ([6]) of multiplicity p, the order of the stabilizer. Moreover all the 
stabilizers are conjugated, therefore they have the same order. Thus we see that all the roots of 
equation ([6]) have multiplicity p and that equation ([6]) is the p-th power of a monic polynomial 
of degree n which we will call U. Let us try to take advantage of the polynomial U and of ^> to 
parameterize a root xi of an anharmonic h n . We remind that the sum of the roots of an anharmonic 
h n can be taken to be zero (by using an homography). We have the proposition 

Proposition 7. Consider one of the t as defined by equation ([5]) and rjj the roots of U . Impose 
the condition 

Xi-.= P(i) -^— + Q(i) 



and that 



with 



P(s(t)) = 



s(z) 



t-rji 
ad — be 

{ci + dy 

az + b 



P(t) 



Then xi admits the parameterization 



ad — cb ^= 0. 



U'(t) 



i - rjj nU(t) 
*'(t) 
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Proof. As 



one has 



But 



hence 



X] Xi = °> 

i 

i \ i 

V 1 _U'(i). 

Q® = 



= 0. 



C/(t)' 
U'(t) 



nU(t) ' 

Now if P and P\ are two rational functions satisfying the conditions of the proposition, then 

Pi 
P 

is left invariant by all linear fractional transformations, therefore it is a rational function with 

p 

numerical coefficients of ^(t). As the latter is supposed to be equal to T e F, — e F. Pi differs 

from P only by a multiplicative factor belonging to F. Suppressing this factor is equivalent by the 
below representation of a;,: 



Xi := P(i) 



1 



t-Vi 



Q(t) 



to multiply Xi by an element of F. This does not change any of the fundamental properties of the 
anharmonics (the Constance of the cross-ratio and the irreducibility) . Therefore it suffices to get 
an arbitrary P(t) in order to solve the problem. Let W be the derivative of the absolute invariant 
with respect to t. As 

'ai+b s 



then 



Now we just take 



at + b 

ci + d 



,ct + d 
*'(t) 



*(t), 



di 



d 



P(t) 



at 



ct + d 
1 



= *'(t) 



(ct+d) 2 
ad — be 



*'(t) ' 



This completes the proof. 

One has therefore for expression of the Xi 



□ 



(8) 
with 



X ■] 



U'{t) 



Vi nU(t) 



*'(t) 



f(t,%) 



*(t) = T and Ufa) = 0. 

We show that the general solution of the Riccati equation admits a representation in the form 
u = f(t, C), C an arbitrary constant. 
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Lemma 8. The general solution u of the given Riccati equation ([3|) can be written in the form 

u = f(i,C). 

Proof. Take three solutions x\, x 2 , x 3 with respective representations 

Xi = f(t,J?i), x 2 =f(t,Tte), x 3 = f(t, 773)- 

Define £: 



u = f(t,p). 

One knows that the cross-ratio of it, Xi, X2, x 3 , is a constant X c 

(u - xi)(x 3 - ari) 



So 



(x 3 - x 2 )(u - x 2 ) 



But by very construction 



(u - xx)(x 3 - x 2 ) (p- - m) 



(x 3 - xi)(u - x 2 ) (m - r}i)(p- %) 
So p is a constant whose value depends on K c . It is the arbitrary constant C. 



□ 



Now we have that 



and 



This leads to 



*(t) = T 
u = f(t,C). 



C = t 



U'(t) 



-1 -1 



nU(t) 



Taking the derivative — , we get 
at 



(9) 



W — h tlW 

at 



d 






U'\ 




{nil) ' 


\- (u&- 




Jt 







= 0. 



This gives the form of the Riccati equation @ associated to the algebraic equation ([3]) all of whose 
roots are solutions of the Riccati equation ©. 



Remark 9. Elimination of t between 

T = *(t) and 



f(t,m 



and setting the resulting equation equal to zero leads to the desired form of the anharmonic equation 
of the theorem ([T]). More precisely, the resulting equation is of degree N; also one sees that the 
construction does not depend of r\i . Moreover the stabilizer Si of r\i (which is the stabilizer of Xi), 
fixes globally the two equations; therefore the equation we get, is the p-th power of a polynomial 
F(xi,T) (as any of its roots will have multiplicity p), the polynomial which we are looking for 
(after an eventual division by the leading coefficient). This monic polynomial of degree n is in 
¥[x] and is irreducible as its p-th power does and its n roots are given by construction by the n 
solutions of the anharmonic h n . 



10 AHMED SEBBAR AND OUMAR WONE 

2.1. Degree of the anharmonics. We have seen that each anharmonic h n , has a Galois group 
H of cardinal TV = np, with p the cardinal of a stabilizer of a point (the stabilizers are conjugate) 
and p\(n — 1) or p\(n — 2). Moreover H is either cyclic, or dihedral, tetrahedral or octahedral or 
icosahedral. Also given any such group H, all its subgroups under the previous constraints are 
susceptible of being a stabilizer of a root of an anharmonic, as in the reasoning for the construction 
of an anharmonic, the latter were arbitrary. 

We are going to use these specifications in order to determine the possible degrees n. 

(1) First of all one remarks that if p = 1, H can be any of the finite subgroups of PSL(2,C) 
and n assumes an arbitrary value > 4, for the cyclic case; an arbitrary even value > 4 for 
the dihedral case. It takes the value n = 12 for the tetrahedral case, the value n = 24 for 
the octahedral case and the value n — 60 for the icosahedral case. 

(2) Next assume p > 1. 

(a) Then H can not be cyclic. This is because the group G acts transitively on the roots 
of h n , so the stabilizers of the roots are conjugate. The group being cyclic, hence 
commutative, we see that all the stabilizers are cyclic subgroups and are the same: 
H Xi = H Xj for i 7^ j. Therefore H Xi fixes all the roots of the anharmonic h n , n > 4 
and is different from the identity (p > 1); this is absurd because of the Constance of 
the cross-ratio (any subgroup of H which fixes more than two roots, fixes all of them 
and H Xi should be the identity). 

(b) If H is dihedral then p = 2; the case p\ (n — 1) gives an odd n and p\ (n — 2) an even n. 
To justify this assertion we remark that the dihedral group group of order N — 2m is 
generated by two transformations 0, e with 

9 m = Id; e 2 =Id, e- 1 9e = e- 1 . 

It contains the 2m transformations 

6\ i € {0, • • • , m - 1} and e9 r , r e {0, • • • , m - 1}. 

We claim that the stabilizer of a root Xi, H Xi can not be one of 0\ i € {0, • • • , m— 1}. 
The argument is based on the transitivity of the Galois group H, whigh gives that 
the stabilizers H Xi are all conjugated. 

If H Xk is of the form (9 ( ) and is conjugated to H Xj with I ^ j, via some eO d , then a 
quick computation gives 

&- d e - 1 e p ee d = e- d Q- 1 e d = q- 1 . 

Therefore H Xj = (6 _1 ); thus because the operation of conjugation is an automor- 
phism, one has §(H Xk ) = §(H Xj ) = j}((6)) = p = m. Hence because N — np = 2m, 
one gets n = 2; absurd. So whenever H Xk = (0 P ) has a conjugated stabilizer via an 
element of the form e® d , we are led to a contradiction. 

We must under this fact examine the other alternative which is: the remaining stabi- 
lizers are all conjugated to H Xp via elements of the form Q d . The conclusion in this 
latter case results from the arguments given when H was cyclic. To summarize, in 
the dihedral case, the stabilizers are of the form H Xi = (eQ l ) and p = 2. 

(c) In the tetrahedral case N = 12 = np. The tetrahedral group contains only elements 
of order 2 or 3. If p = 2, then n — 6. The only case which can happen also in the 
tetrahedral case is p — 3 and n = 4. 

(d) If the Galois group H is octahedral, then N = 24 = np. The octahedral group contains 
elements of order at most 4. lip = 2, then n — 12; p = 3 gives n — 8 and finally p = 4 
corresponds to the value n = 6. 
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(e) Lastly the icosahedral has three cases; p — 2 and n = 30. There appears also the case 
p = 3, n = 20 and finally the case n = 5, p = 12. The icosahedral group contains as 
is well-known only substitution of order 2, 3 and 5. 

These various cases according to p = 1 or p > 1 are summarized in the following tables 



Group 


n 


P 


Cyclic 


n 


1 


Dihedral 


n 


1 


Tetrahedral 


12 


1 


Octahedral 


24 


1 


Icosahedral 


60 


1 



Table 1. The case p = 1. 



Group 


n 


P 








Dihedral 


n 


2 


Tetrahedral 


{4,6} 


{3,2} 


Octahedral 


{2,3,4} 


{12,8,6} 


Icosahedral 


{2,3,5} 


{30,20,12} 



Table 2. The case p > 1. 



Example 10. Assume we are in the cyclic case, then p = 1, n = N, ^ cyc = t™ = T. Also 
[%) =rf-K, with if e C. One has 

x = f(t, V )= ' ! 

this implies 



ni n-i 77 t n -/r 



This can be again written 



r l = l-(nTx + ^- l ) TeW 



i =x 



as such a transformation preserves the fundamental properties of anharmonics. In conclusion we 
get 

T 

3. RlCCATI EQUATION, DARBOUX-HALPHEN SYSTEM AND CHAZY EQUATION 

Our goal in this section is to show that Riccati equations are not exceptional in their property 
of satisfying algebraic equations. Indeed we are going to exhibit a third order polynomial, all of 
whose three roots satisfy a system of Darboux-Halphen type, which in the generic situation is 
irreducible (when the three solutions of the Darboux system are distinct). We will also consider a 
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remarkable irreducible third order polynomial whose three roots satisfy the Riccati equation. We 
fix once and for all a generic lattice in C. 

Let us introduce the basic definitions about modular forms. We recall there is a natural way 
to identify modular forms with some functions defined on lattices (two dimensional Z-submodules 
of C) in C having some invariance and homogeneity. We refer to [2(| for details. Let 9t be 
the set of lattices in the M-vector space C and 971 = {(wi,^) G C* 2 , r = > 0}. Let also 
$j = {z £ C,3(z) > 0}. Then 9* identifies to the quotient £0t/SX(2, Z). Moreover C* acts on £fl 
and on 971 with two more identifications 

m/c*^s), <n/c* «£/psx(2,z). 

A function F : 9\ i-> C is of weight n if F(XA) = \-' l F(A) for all lattices A and all A E C*. For 
(wi, L03) G 971 and A = A(wi, w 3 ) = Zwi © ZW3, the lattice generated by wi and W3, we write simply 
F(uj%, U3) = F(A) so that F(Xuji, Xuj^) = \~ n F(u)i, W3). Moreover W3) is invariant under the 

action of SX(2,Z) on 971 (as the lattice A(wi,W3) is). This implies that there exists a function 

f:Sj^C, F( Wl ,W3)=w 3 -"/( — 

V^3 

far + (3\ 

The invariance of F under SX(2,Z) means that fir) = vyr + S) n f\ r- . For even n, we 

\ 7T + J 

recover the classical invariance property of elliptic modular forms. Moreover as I2 and — 12 act in 
the same way, one can consider that it is the modular group PSL(2, Z) which acts on the functions 
F. 

On the other hand, one of the most general principle in elliptic function theory is that for any 
elliptic functions <^>(m, wi, W3) of periods 2^1, 2u;3, the two functions 

dcj) d(f> d<f> 
f{u) = CJi- h lj 3 - h u— 

OLOl OLO3 ou 

d<p d(j> d<j) 

g{u) = ??i ~ — \-m^ — Kw^- 

ou>i 00J3 ou 

are also elliptic with the same periods 2wi,2oj3. We set W2 = ui\ + W3 and adopt the classical 
notations for Weierstrass elliptic functions 

(10) P(u) = p(u;(Ji,u 3 ) = \ + V / 

m 2 +n 2 7^0 

Let 



(m + 2mwi + 2naj3) 2 (2rawi + 2nw3) 2 



C(ti,wi,w 3 ) : = i + V (- — — 5—— - + 



u V (u + 2mu)i + 2nw3) (2rawi + 2nw3) (2mwi + 2na;3) 2 

rn 2 +n 2 y^0 

One has 

(11) C' = -p, % = C(^i) 5 % = C(w 3 )- 



We set 



and 



(12) 



— ,3r>0, q = e 2i7rT 



02 = 60 ~ — 7. 

,g 3 = 140 Yl ' 



(2mwi + 2nw 3 ) 6 ' 

m 2 +n 2 =£0 V 1 A > 
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As a result 

dp dp dp 

-2p=u x - hw 3 7 — ^ U ^T 

, s dui dw 3 du 

„ 2 1 d<t> d<f> ^d<t> 

- 2p + s 92=m d^ + m d^ + (:{u) d^- 

The constants g2, g% are related to the classical Eisenstein series £4 and E§ by 

4 °° „3„n 



(14) 



Set also 



We define 



E 4 (r) =12(^1) g 2 = 1 + 2405: 

6 00 

£ 6 (r) =216 (^) 53 = l- 504^^ 



(7" 

n=l y 

OO ^ n 

rrq 



1 - 24 5] cri(n)g" = 1 - 24g - 72g 2 - 96g 3 - 168g 4 



Definition 11. Let T be a subgroup of finite index of the the modular group SL(2^1^). A mero- 
morphic (respectively holomorphic) function f : fj h->- C is a weak meromorphic modular function 
of weight k with respect to T if 

for every \ . ) € <SX(2,Z). // / is meromorphic (respectively holomorphic) at the cusps of T , 



7 ^ / 

zi wiZZ &e called a meromorphic (respectively holomorphic) modular form of weight n with respect 

to r. 

These above Eisenstein series are modular forms of weight 4 and 6 respectively for the the 
subgroup SL{2, Z). It is a fundamental result, that will be used later, that every modular form for 
SL(2, Z) is uniquely expressible as a polynomial in £4 and Eq and the extension ring C[£2, £4, Eq] 
of C[i?4, Eq] is a differential ring. More precisely, the following basic relations of Ramanujan hold 

(15) ^-T"^6 = \{E 2 E & -Et) 

Ztir dr I 

1 d 1 2 

2^d^ 2 = 12^-^). 

In other words, the field ( C (E2, E4, E§) , — ) is a differential field. The subfield of constants 

\ 2l7T czt J 

is C (as it embedded in the field of meromorphic functions on fj). 

For the particular value u = Ui, i = 1, 2, 3, the equations (| 13[) become 



(16) 



det dei 

-2et = uji- 1- uj 3 - — , 

dui\ dojz 

1 dei dei 



duii du>3 
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We will later deal with the two partial differential operators 

D = — + — 

2 771 doji ^ 3 dais 

Its importance lies in the fact that it converts ellipticity properties into differential relations for 
certain modular forms HI, HI, U|. 

3.1. ei, e 2 , e% and Darboux-Halphen system. 

Theorem 12. For every r, St > 0, the function e,-(r) = p{uii; 2u>i, 2uj 3 ), i = 1,2,3, solves the 
nonlinear differential equation of Riccati type, with coefficients in C (i? 2 , E4, Eq) 

V = % - (-V 2 + y ^ 2 (r)y + ^ 4 £ 4 (r) 
Proof. We have seen in equation (fTBj) . that for the operator D 2 = —2r\\-J^ — 2r\ 3 -J^- 

D 2ei = 4e? - |g 2 , i = 1,2,3. 

Ul 3 

We consider the new independent variables u>i, r = — . The identities 

e*(2wi,2w3) = (2wi)- 2 e fc (l,r), fc (2wi,2w 3 ) = (2wi)- 4 <? fe (l, r), 

7T 2 c / s 71-2 2 7-, / \ inT 

VlUJl = 12 7,2 W2 = 12 T 2 ^ ~ "T" 

and a straightforward calculation give the theorem. □ 

The partial differential operator D2 can be replaced by a single differential operator. For u) = 
2iTTT, ei verifies 

(17) A (w 2 ei ) - Q. 92 w 4 + ^—e^\E2 ~ 2e 2 w 4 ^ , 1 < i < 3. 

The system of differential equations (|15p is now necessary to show that the system (|17[) is leads to 

4i 

a Darboux-Halphen system. Setting t = — r and 

7T 



the system ([T5|) becomes 



.T = £Jo, W = £m, Z = Ek, 

12 ' y 12 ' 216 ' 



dx 1 2 1 



(18) %=2xy-3z 

dt 

dz 1 2 

dt 6 y 

With Xfc = Y2-B2 + jGk,k = 1,2,3, the equations (1171) take the form 

j t {X 1 +X 2 )=X 1 X 2 

(19) l(x 2 +X 3 ) =X 2 X 3 

dt 

^-(X 3 +X 1 )=X 3 X 1 , 
dt 
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which is a Darboux-Halphen system. Moreover It is transformable into (|18l) by means of the 
substitutions: 

x=~(X 1 +X 2 +X 3 ) 

(20) y = t (Xf + xl + X\ - X X X 2 - X 2 X 3 - X 3 X X ) 

z = A (2X a - X 2 - X 3 ) (2X 2 - X 3 - Xx) (2X 3 - X! - X 2 ) . 

3.2. Chazy Equations. A particular case of Chazy equation is the following third order differ- 
ential equation 

(21) 7"' = 677" - 9 7 ' 2 , 7 = 7 (r). 

It has the particular solution 7(r), with j(t) = X^n>o a nQ n > 1 = e 2i7rr , unique up to the substitution 
TH>T + T .a„i-} a n e 2mTa . It is related to the Eisenstein E 2 as follows 

7(r) - - 24q - 72q 2 - 96q 3 - lb8q 4 -•••) = -\e 2 (t). 

b b 

This can be seen from (TT5t by eliminating E4 and -E^- It is a remarkable fact [nj that if 7 solves 
the Chazy equation, the three solutions wi(r), uj 2 (t). uj 3 (t) of the cubic equation 

(22) ^ + ^{t)^ + ~ 7 >(t)u> + iy'(r) = 
are solutions of the Darboux-Halphen system: 

LJl = -U>l{uj 2 + UJ 3 ) + UJ 2 UJ 3 

(23) W2 = -w 2 (wi + w 3 ) + 

w 3 = — w 3 (w\ + w 2 ) + wiw 2 . 

which is equivalent to (fT9|) by setting X, = — 2wj, i = 1, 2, 3. 
Proof. One has 

wi +cj 2 + w 3 = --7(r), 
3 

and 

cjicj 2 o; 3 = -^7"(t). 

Derivation of the three relations with respect to r gives 

3 

lo x = -u> 2 -cb 3 - -^'(r), 

3 \ 3 
-lu 2 - uj 3 - -7'(r) (uj 2 + lj 3 ) + cj 2 wi + cj 3 wi + cj 2 w 3 + lu 3 lo 2 = -7"(r) 



and 



3 \ 1 

-ci>2 - W3 - 2^'^ ) W2 ^ 3 + ^2^i^ 3 + w 3 u;iu; 2 = — -7"'(r). 
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This gives the system 



• ■ 3 , 

UJl = -uj 2 - 0J3 - -7 



3 3 

wa(wi - lu 2 ) + tl) 3 (ui - W3) = -jl" + -jl'{u 2 + oj 3 ) 

1 3 

oj 2 uj 3 (oji - uj 2 ) + uj^uj^ujx - uj 3 ) = -jl'" + -jl'u 2 uj 3 

Set A = (wi — u 2 )(uji — W3)(w2 — W3); then suppose the three roots are all distinct; applying the 
Cramer formula one gets for the value of ui 2 , the following 

13 3 

jY'i^i - W3) + ^'W^i - cj 3 ) + -7'w|(wi - cj 3 ) 

w 2 = £ • 

Using the fact that 7 satisfies the Chazy equation 

7"' = 677" — 9"/ 2 , 7 = 7(7-) 

and the expressions of 7' and 7" in terms of wi, CJ2 and CJ3, one gets after simplification the value 
of Cj 2 

oj 2 = —u 2 (lui + 0J 3 ) + (J1UJ3. 
The values of lo\ and L03 follow in the same manner. □ 

Remark 13. Let Y — PSL(2,1) be the modular group and Y(2) be its principal congruence 
group of level 2. The corresponding Riemann surfaces Sj* /Y = (S) U {oo})/T and Sj*/T(2) = 
(fjUjoo, 0, l})/r, after suitable compactification, (adding cusps) have genus and (after a suitable 
normalization) elliptic fixed points of orders 2 and 3 at i and e %% ^~ for the case of Y . Denote by 
J and A holomorphic mappings which establish isomorphisms H* /Y = P 1 and H* /Y(2) = P 1 . 
These "Hauptmodules" in the terminology of Klein, are modular functions for Y and T(2). One 
can express the general solutions of the Darboux-Halphen system and the Chazy equation in terms 
of these functions. We have [2(| 21 , 28 1 that the transformation rule 

>/ n 1 f ar + b\ c ( a b \ „„,/_ ^ 

w ( T ) -> 7 i^wi : + ;> 7 )ePSL(2,C) 

v ; (cr + d) 2 \cT + dJ cr + d \ c d J y ' 



together with the formulas 



Id, A' 
Ul = 2 dr T 
Id, A' 



2 dr A — 1 
Id, A' 
W3 = -o"7~ lo g" 



2dr °A(A-1) 

provide the general solution to the Darboux-Halphen system. For the case of the Chazy equation 
we have that the transformation rule 

1 (aT + b\ 2c ( a b\ nOT ^ 

and the formula 

1, (J') 6 
7= 6 l0g JV-l) 3 
provide the general solution for the Chazy equation. 
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When the three solutions of the Darboux-Halphen system are not distinct, for example when 
CJ2 = W3, say, then 

c a c 



ct + a (cr + ay 



When they are all equal they are given by 

1 

U)l = L0 2 = 0J3 = 



with To an arbitrary constant and modulo the transformation property. 
The Chazy equation admits also the rational solution 

2c 2 a 



7 



cr + d 3 (cr + d) 2 ' 



Remark 14. The Darboux-Halphen system arose historically in 1878 in the study by Darboux of 
the existence in Euclidean space of a one-parameter family of surfaces orthogonal to two arbitrary 
given independent families of parallel surfaces (such a family is necessarily quadratic and ruled), 
and it was solved in 1881 by Halphen. This equation occurs in the Bianchi IX cosmological 
model and also appears as a special reduction (obtained by change of parameters) of the self-dual 
Yang- Mills equation. See (2^ . 

3.3. The degree three polynomials for the Darboux-Halphen system and the Riccati 
equation. We remind that the Weirstrass function 

p(u) = p(u;uj 1: uj 3 ) = \ + V 

m 2 +n 2 ^0 

satisfies the equation 



(u + 2muj\ + 2nuj 3 ) 2 (2mwi + 2nu 3 ) 2 



(24) p' 2 =4p 3 - 52 p-ff3 
But we have 

(25) p <(l)^(I). p .(l±I). 0i 

therefore with the previous notation of theorem (|12[) . we see that e±, e 2l e%, solve the equation 

(26) AX 3 - g 2 X -g 3 = A(X - ei )(X - e 2 )(X - e 3 ) = 0. 

Remark 15. One can show that e±, e 2l e 3 are modular forms of weight 2 for T(2) and that they 
are distinct as the discriminant of the non singular elliptic curve Y 2 = 4X 3 — g 2 X -~ g 3 is the 
discriminant of the degree three polynomial AX 3 — g 2 X — g 3 up to non zero constant factor. 

We have the well-known fact 
Proposition 16. For all ^ ^ ^ ^ G SX(2,Z), we have 

E 2 (^^4) = (7T + S) 2 E 2 (t) + 6^ ( 7 r + 5) . 
\yr + o J itt 

Now as ei, e 2 and e 3 satisfy the identity 

ei + e 2 + e 3 = 0, 
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one easily sees that the e% can not all belong to C(E 2 , Ea,E§) as the functions E 2 , E4 and Eq are 
algebraically independent over C. Our goal is to show that the polynomial 

x*- 9 -lx- 9 -l 
4 4 

whose three roots are the € M2(T(2)) is irreducible on the differential field C(E2, -E4, E§). g 2 
is modular of weight 4 and g% modular of weight 6 for SL(2, Z). Now we know that the are 
modular only of weight 2 for T(2); to show that 

X 3 - 9lx - — 
4 4 

is irreducible on the field C(E 2 ,E4, Eg), it suffices to show that it is integral on C[E 2 , E4, Eg] and 
then use Gauss irreducibility criterion. This is equivalent to the fact that none of the belongs 
to the latter field. To see this is true, we assume for instance it were the case that e\ belonged 
to C[E 2 , E4, Eq]. e\ is modular of weight 2 for T(2), using the transformation rule for E 2 applied 

to the matrix .9=^2 l J (proposition [T6|) . one sees that E2 is not modular for T(2). As 

the weights of E4 and Eq are 4 and 6, e\ can not be a polynomial in them. Thus one gets a 
contradiction. Therefore we have the proposition 

Proposition 17. The solutions ej of the Riccati equation of theorem are the roots of the 

irreducible degree three equation X 3 — — ^ on the differential field {E 2 , E4, Eg) , 

The Galois group of the equation is thus S3 . 

We want now to show that in the generic case (all three solutions distinct) the solutions 

of the Darboux-Halphen system are algebraic over the differential field (0(7,7', 7"), — | C 

V dT J 

C(E 2l E4, Eq), ~J~^J > w ith the solution of the Chazy equation also considered in the generic case. 
For that we introduce the concept of quasimodular forms due to Zagier and Kaneko [l6j . 

Definition 18. A (meromorphic) quasimodular form of weight k and depth p onT is a meromor- 
phic function f on Sj such that 



(27) 

_ / a b , 
zef),{ _ d ] er 



c 



and where the fi are meromorphic on f) with moderate growth at the cusps. The space of quasi- 
modular forms of weight k and depth p on Y is denoted by M { k ~ p) = Nlf p) {Y). The prototype of 
quasimodular form is the Eisenstein series E 2 that is of weight 2 and depth 1. 

The following summarizes the properties of quasimodular forms 

Proposition 19 ([16]). Let T be a modular subgroup and let k and p be non negative integers 
• The space of quasimodular forms on T is closed under differentiation 
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• Every quasimodular form on T is a polynomial in E2 with modular forms as coefficients. 
More precisely we have 

Mp) = (B p r=Q M k „ 2r (r).E r 2 

for all k, p > 0, where Mj(T) denotes the space of weight j modular forms on T. 

• Every quasimodular form on F can be written uniquely as a linear combination of deriva- 
tives of modular forms and of E 2 . More precisely we have 

Mi<P) = h P r=0 D r (M k - 2r (T)) ifp < k/2, 

k \® k r =o~ lDr ( M k-2r(r)) ® £D k ' 2 - 1 E 2 ifp > k/2. 

X' X' X' 

In the generic situation, the functions f\ = — , f 2 = and fa = — — (remark [13]) are 

A A — 1 A(A — 1) 

modular of weight 2 for T(2). Therefore their logarithmic derivatives, hence the u>i are quasimodular 

forms of weight 2 and depth 1 for T(2) (just take the derivative of the modular relation). Moreover 

1 A' 

E 2 = ^^"^"j with A = q Yl n>1 (l — q n ) 2i the modular discriminant of weight 12. Thus taking into 
account the previous theorem (lq |. one sees that none of -r can belong to CE 2 as this would lead 

Ji 

to a relation 

fi A 

with c e C*. This is a contradiction as /j is of weight 2 and A of weight 12 for T(2). In other 
words one has a non trivial relation 

f. 

-j- = 9i + CiE 2 

fi 

with gi £ M 2 (T(2)) and Cj ^ 0. One immediately sees, using a similar approach to the case of 
/' 

proposition ITB|) . that -j- ^ C[E 2 , E4, Eq] as this latter contains no modular form of weight 2 for 

fi 

T(2). We reach therefore a similar conclusion to the case of proposition [16] that is 

Proposition 20. If the solutions of the Darboux-Halphen system are generic (ie all distinct) and 
also if the solution of the Chazy equation is generic (not a rational function) then the polynomial 

uj 3 + ^(tP 2 + |V(t)w + \i\r) 

is irreducible on the differential field ( 0(7,7', 7"), — J • Therefore the extension C(w)/C(7, 7', 7") 
is galoisian with Galois group S3. 

4. Algebraic solutions of the Riccati equation on C(p(z),p'(z)): Another approach 
In this section we investigate the following question: for which potential q the Riccati equation 

u + u = q 

admit an algebraic solution over the field C(p(z), p'(z)). We will discover this is given by the 
anharmonic Weirstrass function po corresponding to the case g 2 — 0. 

First of all we recall some notions from classical invariant theory. Let if be a field of character- 
istic zero. A binary form is by definition a homogeneous polynomial 

f(x, y) = a x n + r)aix n ~ 1 y+ r\a 2 x n - 2 y 2 + ■ ■ ■ + a n y n 
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with coefficients in K, n is the degree of the form /. Another used notation is 

f(x,y) = (a ,a 1 ,---a n )(x,y) n . 

The set V n of binary forms of degree n is a if -vector space of dimension n and can be identified 
with the space of polynomials K[a , a\, ■ ■ ■ a n ] on which the groups GL(2, K) acts in the following 
way: 

a b \ 

eGL(2,K),feV n , (gf)(x,y) = f(ax + by,cx + dy). 

c d J 

We introduce the following important differential operator 

n 92 92 



dxdy' dx'dy 

It is known as Cayley Omega process. 

For two given binary forms Q(x, y), R(x, y), their transvectant of degree r is the function 



t d r Q d r R 



d r ~ l xd l v d l xd r ~ l v 

The transvectant of degree r can be obtained from the Cayley Omega process using the following 
formula 

(Q,R) r = n r {Q(x,y),R(x',y')} lx ,= x , y ,= y . 

For example 

(Q > R 1 ) — QxRy QyRxi (Q i R*) — QxxRyy 2Q X yR X y -\- QyyR X x~ 

X 

For polynomials of one variable F, we consider the projective coordinate p = — and define the new 

V 

x 

polynomial Q(x, y) = y n F(—). It follows that if F, G are polynomials of degrees n, m respectively, 

V 

the r th transvectant of F, G is, for r < min(m, ri) 

k (r\ (m - fc)! (n-r + fc)! d T '- k F d k G 



k J (m — r) ! (n — r) ! dp r k dp k 



(F, G y = J2(-i 

It is of degree < m + n — 2r. Few examples are given by 
(F,G)° - FG, 
[F,G) X = mFpG-nFGp, 

(F, G) 2 = m(m~l)F pp G-2(m~l)(n-l)F p G p +n(n-l)FG pp . 
and the Hessian 

H(F) = \ (F, Gf = n(n 1) (fF pp - . 

We need in the sequel the fourth transvectant 

(F, F) A = 2(m- 3)(m- 2) (m(m - 1)F (4) F - 4(m - 3) (to - 1)F (3) F' + 3 (to - 3) (to - 2)F" 2 ) 

In particular the fourth order differential equation corresponding to the vanishing of the fourth 
transvectant is 

(28) to(to - 1)F (4) F - 4(to - 3)(to - 1)F (3) F' + 3(m - 3)(m - 2)F" 2 = 0. 
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Following [22] the equation (|28|) can be reduced to generalized Chazy equation which is a third 

........ F' 

order differential equation. For if R is given by —mR = — , then an easy calculation gives 

F' = -nRF 

F" = {-nR' + n 2 R 2 ) F 
(29) = (-nR" + Zn 2 RR' - n 3 R 3 ) F 

= (-nR® + 3n 2 R' 2 + in 2 RR" - 6n 3 R'R 2 + n 4 i? 4 ) F. 
So that (F, F) 4 — becomes the generalized Chazy equation for R 

R^ - 12RR' + 18i?' 2 = (R' -R 2 ) 2 . 

n — 1 

On the other hand the coefficients of the fourth transvectant (/, /) can be computed recursively 
from those of the binary form /. For f(x, y) = (do, ffli, a 2 , • • • , ci n )(x, y) n we have 

\(f,f) i (^,y) = (ao,ai,a 2 , • • • ,a m )(x,y) m , m < 2(n - 4) 



with 



(30) 



ceo = aoa4 — 4ai<23 + 3a 2 

m\a.\ = (n — 4)(aoa5 — Zaia^ + 20203) 

(n — 4)(n— 5). 2 . , . N o, _ 2n 

m 2 a 2 = (aoa6 - 4aia 5 + 7a 2 a 4 - 4a 3 ) + [n — 4) (aia 5 - 4a 2 a 3 + 2a 3 ) 



a m — a„a„_4 — 4a„_ia n _3 + 3a^_ 2 



"Jm-iOm-i = (n - 4)(a„a„_ 5 - 3a„_ia„_4 + 2a„_ 2 a„_ 3 ). 
More generally according to r is even or odd: 



m. 



~l r QL r — ^ ^ Pr,sFr,s: Tfl T OL T — ^ ^ Pr.sFr. 


where 

(31) m r = ( m ), p r = ( n 4 ) ( n 4 

\r ) \ s J \r — s 

and 

(32) P rs = a s a r - s +4 — 4a s +ia r _ s+ 3 + 6a s+2 a r _ s + 2 — ^a-s+3 a r-s+i + a s +40r-s- 

We now define on the ring C(p(z), p'(z))[u], the following differential operator, which will play a 
important role in the sequel 

X:C(p(z),p'(z))M — t C(p(z),p'(z))M 

Let <I>(u) € C(p(z), be the irreducible polynomial of minimal degree n satisfied by a 

solution of the Riccati equation: vl + u 2 = q. The assumption of minimality means the irreducible 
polynomial of least degree satisfied by an eventual solution of the Riccati equation. 

(34) $(u) =u n + Y\ aiUn ^ + n(n 2 7 1) Q2 ""' 2 + ■ ' ' + a « 
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then 



(35) 



X($) 



n(ai — m)$ 



We define 



Definition 21. Let K be a field of characteristic zero, K[X\ ■ ■ ■ ,X n ] a ring of polynomials over 
it and D a derivation on K[X\, ■ ■ ■ ,X n ]. A non-zero element P of K[X±, ■ ■ ■ ,X n ] is a called a 
Darboux polynomial for D if P is an eigenvector of D for its action on K[X\, ■ ■ ■ ,X n ], namely 
there exists R £ K\X\, • • • , X n ] : 



When R = P is called a first integral. The eigenvalue is called a Darboux cof actor. Two Darboux 
polynomials having the same cofactor give rise to a first integral. 

So <5>(u) is a Darboux polynomial for the mentioned derivation, ie a non trivial polynomial 
eigenvector of the derivation 



Remark 22. It is perhaps worthwhile to observe that in general, an algebraic function f on a 
differential field (L,5) of characteristic zero verifies a linear differential equation with coefficients 
in L. In fact if /" + <2i/ n_1 + • • • + a n = 0, a; 6 L, 1 < i < n is the minimal polynomial of f over L, 
then necessarily n/"- 1 + (n- l)ai/™~ 2 + • • • +a„_i ^ and /' e L(f). Hence /", / (3) , - •• € L(f). 
The dimension of the the L-vector space L(f) is n and therefore the vectors f^ n \ / < -™~ 1 - ) , • • • , / of 
L(f) are linearly dependent over the field L, that is to say f is a solution of an homogeneous linear 
differential equation. 

Remark 23. As our base field K = C(p(z), p'(z)) is a field of meromorphic functions, the existence 
of a non-trivial first integral for X, P(u), is equivalent to P(u) being constant on the solutions. 
Therefore the general solution of our Riccati equation is given in this case by a relation 



£ an arbitrary constant belonging to C. Moreover two first integrals for X can not be functionally 
independent as we know from the Cauchy theorem for differential equations, that the general 
solution of a first order equation depend only on one arbitrary parameter. 

We have the following well-known lemma 

Lemma 24 ([8]). A non trivial polynomial P(u) in C(p(z), p'(z))[u] is a Darboux polynomial for 
d d 

the derivation — h (q — u 2 )— if and only if all its roots are solutions of the Riccati equation 

oz du 

du „ 

— + u z = q. 

dz 

Using equation (|35|) . one sees that the coefficients of verify the following relations 



D(P) = RP. 



v 9 , 2n 9 
X= dz + {q ' u) d-u 



P(u) = € 



(36) 



(n — l)a2 = na 2 ~~ a 'i ~ Q 

(n — 2)d3 = na\a2 — a' 2 — 2a\q 

(n — 3)a4 = na\a^ — a'^ — 3a 2 q 



(n - k)au+i = na x a k - a' k - ka k -iq. 
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This shows at least that if it happens that a\ , q belong to some differential subring L C K then so 

do the other coefficients a,k- 

We consider the homogeneous polynomial 

K( Ul ,u 2 ) = — ) 

U2 



UA. 



k=0 



and its Hessian 



2 K 
du\ 

8 2 K 
du\du2 



d 2 K 
du\du2 

8 2 K 



= (n - \)u 2 2 n - A (n$$" - (n - 1 )$' 2 ) . 



The polynomial H — n<I><I>" - (n — 1)$' 2 ('= J^) is of degree < 2(n - 2). For latter use we have 
the following important property 



Lemma 25. Let K be field of characteristic zero and = u n 
... + a n an irreducible polynomial then H($) ^ 0. 

Moreover the partial differential operator X is a derivation 

X(uv) = uX(v) + vX(u) 
with the following commutation property 



ai« + — — L a 2 u n - 2 + 



I! 



2! 



Hence 



and 
(37) 

The Jacobian is 



d d 



X(&) = -n$ + (nai -m + 2u)$' 

= -2(n - f )*' + (nai + 4u)$" 
= -2(n - 1)$$' + (nai -™ + 2u)$$" 
X($' 2 ) = 2$' (-n* + (nai -™ + 2u)*') . 

X(i2") = 2(noi - riw + 2u)H. 



3£ 


d± 


dx 


du 


OH 


dH 


dx 


du 



X(H) H' u 



nai{$H -2H<f> ) - uQ. 



with n = n$i?' - 2(n - 2)H& . It is easily seen that 

X(H') = -2(n-2)H + 2(n ai -nu + 3u)H' 

X{$H') = -2(n-2)$i? + 2(nai - nu + iu)<$>H' + n{a x - u)<$>H' 

X(H<5>') = -n$H + 3(nai -nu + 2u)&H. 
from which we obtain 

(38) X(n) = 3(n ai -nu + 2u)n 
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and finally 

Or 
(39) 



3X{H)n = 2HXVL 



H 3 



= r. 



where T is an arbitrary constant in C (see remark | 
constant. We consider the Jacobian of $ and ft 



Now we are going to exhibit another 



0® 


o± 


dx 


du 




m 


dx 


du 



X(4>) K 

x(n) n' u 



= rwii($n / -3n$') - wOi. 



with X = n$0' - 3(n - 2)0$' 



In a similar way as before, we have 

X(fl') = -3(n - 2)0 + (3nai - 3nu + 8u)Q' 

= -3(n - 2)$0 + 4(nai + 2u)$0' 
X(Q$') = -n$0 + 4(nai - nu + 2u)$' 

so that 

(40) X(Cli) =4(nai -nn + 2u)n x . 

The two identities ([3"7]). (gDJ give 

20 1 X(iJ') = X(fii)F 

that is, with another constant a 

As we said before in remark [23] two non trivial first integrals are functionally dependent and this 
should force differential relation for our $(u). Indeed one has the following 



H = n$$" - (n - 1)$' 2 

= - (n - 2)$'$" 

Q = n$iJ' - 2(n - 2)#$' 

= n 2 $ 2 $'" - 3n(n - 2)$$'$" + 2(n - l)(n - 2)$' 3 

O' = -n(n-6)$$'$"' + 3(n-2) 2 $' 2 $"-3n(n-2)$$" 2 + n 2 $ 2 $"" 

Q x = -4n 2 (n - 3)$ 2 $V" - 3n 2 (n - 2)$ 2 $" 2 + n 3 $ 3 $"" + 12n(n - 2) 2 $$' 2 $" 

- 6(n-I)(n-2) 2 $' 4 

aH 2 = a (n 2 $ 2 $" 2 - 2n(n - I)$$'V' + (n - 1) 2 $' 4 ) . 

In conclusion splitting the two members of the equality Qi = (1H2 into monomials involving only 
powers of $' on one side and the rest on the other side, one sees that $ should divide $' unless 

(n - 2) 2 

a(n — l) 2 = — 6(n — l)(n — 2) 2 or a — —6 . Moreover we have the divisibility of fix — a -ff 2 
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by $ 2 giving the vanishing of the fourth transvectant 

^ 1 = — 

' 4I) = n(n - - 4(n - l)(n - 3)$'$"' + 3(n - 2)(n - 3)$" 2 

= 0. 

a o 2 ^ 2 

Set L(u) = — — — — — . It is a general fact that because — 7 is a first integral of the Riccati 
ou H i 

d d 

equation and \X, — = 2u— that the following relation holds 

ou ou 

(42) = -2uL(u). 

So = L(u)H 6 is a Darboux polynomial. Computation gives 

Li(u) = n(2fm' - 3nH')H 2 ■= loEh 2 . 

5 = 2HQ' — 3QH'. We want to establish a relation between S and $ which will enable us to give 
another first integral of the considered Riccati equation and show that a\ = 0. First we calculate 
S$. One has by very definition 

S$ = 2H<S>n / - 3fl$H'. 

But 

$H' = -(Cl + 2(n- 2)H&) , 
n 

so 

n 2n 2 

S$ = -— + (O + 2(n - 2)ifty) = ft 2 - -i/fii. 

3 3 3 

As fii = aif 2 , with a the previously given value, one gets 

= n 2 - aH\ a = -4- — . 

3 n- 1 

We remark that with this value of a, $ 3 is a factor of f2 2 — aH 3 . 

Let us take the derivation of the previous relation with respect to u; this gives 

-- (H$' + $5') = 2S1^' - 3aH 2 H'. 
3 

One sees that 



2H 



there we obtain the equation 



-- (H$' + $S') = 2 ^ + mH ' - 3aH 2 H' 
3 2H 

= n E + j nH ' -3aH 2 H' 
H 

= ^+ 3 -^(n 2 -aH*) 
= -2(n - 2)S$'. 

Thus 

$s' = 6(n ~ 2) ^s - *'e 

n 
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and this gives finally 

(43) = ( 5 - — ) $'H. 



n 



It follows that 



$£' = ( 5 - — ) S$' := (fc - 1)5$' 



71 

or 

(44) 5 = C^ 1 

where C is an arbitrary elliptic function. 

Now 5 is a Darboux polynomial for X; indeed the previous relations for X(H), X(£Y), X(H'), X(O') 
gives us the relation 

X(E) = (5n(ai-u) + 12u)H 
and this together with the fact that $ is Darboux leads to the following relation involving C 

(45) C = = 12oiC. 

az 

But the primitive of a non vanishing elliptic function for a generic lattice is not elliptic. We refer 
to " 



25| where the various integrals of p n ,n £ Z are computed and related to equivariant functions 



for modular subgroups T C PSL(2,Z). More precisely h : fj — > P 1 meromorphic is equivariant 
under the action of T if 

r+jA _ a/i + 6 / a b \ 
^eFTd) ~ ch + d' \ c d) El - 
This forces cti = 0. The conclusion is with a new constant (3 

n 2 - aH 3 = /3$ fe 

that is to say 

$ fe T-a 



H 3 " f3 

which is another form of the equation p9p . The consequence is that k is an integer satisfying 
k = 6 — — or n(6 — k) = 12. This is because H € C(p(z), p'(z))[u] and 3> is irreducible so any of 
its fractional powers is not a polynomial. We will examine separately the cases n = 2 and n = 3. 
The main result concerning the minimal polynomial $ are given by the following theorem 



,2 



Theorem 26. For § to be a minimal polynomial ()34j) . i£ is necessary that 

a 1 = 0, X($) = -ntt$, r 4 ($) = 0. 

To give an example, we look at the case of n = 4. In this case 

1 1 / 1 /, 2 

3 6 6 

so that the potential q solves an equation similar to (|46p : q" = 

(46) u" + 2cu 2 = 0. 

It has for c ^ 0, the solution 

3 

u(z) = — p (z) 
c 

where po(z) is the equianharmonic Weierstrass function 

z = 
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(n - 2) 2 

For general n = 4, 6, 12 , the potential solves the equation q" = 6aq 2 , with a — — . 

71—1 

From lemma[2Uwe know that a Darboux polynomial for X has necessarily all its roots satisfying 
the associated Riccati equation. As the degree of H is 2(n — 2), we have to examine the cases 
n = 2 and n = 3 separately as in these two cases, the degree of H(u) is presumably less than n. 

• Take the case n = 2. Then 

<&(u) = u 2 + 2aiu + a 2 . 

After computation one gets H(u) = 4a 2 using the condition a\ = 0. But H(u) is a Darboux 

d 

polynomial for X, so (l37l) a, = -^-aa = 0. This forces a2 to be an element of C. As ai = 0, 

oz 

a 2 G C, we get a contradiction as is assumed to be irreducible and C is algebraically 
closed. Therefore the case n = 2 for monic irreducible is impossible. 

• Next consider the case n = 3. Then assumes the form: 

= ?i 3 + 3ai?j 2 + 3a 2 u + 03. 

The expression for H(u) is: 18(a2it 2 + a^u — a 2 ). The previously analyzed case n — 2 shows 
that -ff (w) = 18(a2U 2 + asu — a 2 ) if a 2 7^ necessarily factors on C(p(z), Moreover 
a 2 can not be zero, because otherwise H(u) becomes a polynomial of degree or 1. Since 
we know that is irreducible, we get that is a solution of our Riccati equation (03 can 

dii 

not be as <£>(u) is assumed to be irreducible) which therefore takes the form — — \-u 2 = 0; 

dz 

the general solution of the latter equation is given by u(z) = , Co € C. We claim 

z + C 

that such an u is not algebraic over C(p(z), p'(z)). Indeed if that were the case, then one 
would get a relation of the form 

1 

a 3 = 



(z + Co) 3 

with 03 ^ doubly periodic (minimal polynomial of degree 3). This is clearly absurd. 
Thus a 2 never vanishes. 

For the case n = 3, the identity T4(<I>) = disappears and only the differential system (|36|) . with 
ai = remains 

2a 2 = —q 

(47) (n - 2)a 3 = -a' 2 

= — 3a 2 q 

(3 - 2) 2 1 

which gives q" = 3q 2 which also of the above form with a = — — = - . 

& H y (3-1) 2 

Remark 27. The appearance of the potential q, q" = 6aq 2 or its solution, the anharmonic Weirstrass 
p function is one the unifying point of this study. Despite the similarity, it is very different from 
the Lame potential. 

5. Hypergeometric Functions and Conclusion 
We consider the precedent differential equation for the potential 

(n - 2) 2 

(48) q = 6aq, a = f- 

n — 1 
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with the solution q(z) — — p(z) = — p(z; gi — 0, 173) and we would like to make a connection with 

a a 

some known facts on hypergeometric functions. For commodity reason we set 93 = 4A 6 and the 
Riccati equation for u becomes 

(49) ^ +u 2 = g(z) = I (z) _ 

dz a 

By making the change of variable X — p{z) we obtain 

(50) ^ 2 = 4 X 3 - ff3 = 4(X 3 -A 6 ) 
and if we set 

(51) X = \ 2 ^ Xz = t, u = Xw 
the equation ([4Uf becomes 

(52) -£.«-«.) 
and the equation (|50[) transforms into 



( 53 > UJ =4(e 1} ' 2dt= v^- 

The equation ([52]) is much more simpler than (|4"9")l . 

Now in the equation (|52p we make the change of function 

(54) t = 0(0, V = w 9'-L r 
which transforms a Riccati equation 

(55) ^ +w * = R{t) 
into another Riccati equation 

(56) Tt +v2= r{0 = 6 '" m V 6 ^ 6 ' = I 

so that (|56p is similar to a Darboux-Backliind transformation for Riccati equations. We thus obtain 
a simpler form for the equation (|49)l and (|52|) 



^ , 2 C cqC 3 - 4ci 4 1 

(57) + w = T67F~TF' C0 = a" 3 ' Cl = « + 6 - 

One can show that this last equation reduces to the desired form 

dW , T 2 1 /^A 2 -l !/ 2 -l A 2 -^ 2 + z/ 2 -l 



^ =^ = 4V * 2 ' (I-*) 2 ' .(!-.) 

with 

_ 1 _ n 1 
A ~3' ^-6(^2)' V ~2- 
Before stating the conclusion, we would like to insist on the major role of the hypergeometric 
function in this study. We use classical notations and recall essential facts of the theory of Fuchs 
and Schwarz. The most general form of the hypergeometric equation is 

" (2-A-/i)x + A-l , (1 - A - /i) 2 - f 2 

< 59 > y + ~ — wh) — y + 4,(1-,) y = °- 
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We introduce the new constants a, b, c connected to the local exponents and their inverses by 

a + b + c= l — A — /x, a — b = v, c = 1 — A, 
Mo = 1 - c, hi = c — a — b, fioo = b - a 

fco = , k\ — , /loo = . 
Mo Mi Moo 

The associated Schwarzian equation is 

1 - Mo 1 - Mi 1-/4, 



(60) -4{s,a:} 



x 2 (l — x) x(l — x) 2 x(l — x) 



If (fco,fci,fcoo) is one of the triplets (2, 2,m), (2,3,3), (2,3,4), (2,3,5), the monodromy group of 
the hypergeometric equation is of finite order N = 2m (dihedral), N = 12 (tetrahedral), TV = 24 
(octahedral), N = 60 (icosahedral) , with 

2 1 1 1 

N k fci fcoc 

The main conclusion is that the cases for which the considered Riccati equation has algebraic 
solutions correspond exactly to those cases for which the hypergeometric equations has algebraic 
solutions. These correspond to platonic solids or regular polyhedra. 

This is a general fact, indeed Baldassarri, Dwork and also Maier, following Klein, have shown 
@> S 0> S El that any algebraic solution of a linear second order differential equation with 
regular singularities on an algebraic curve, is the pull-back via a rational map of a particular set 
of hypergeometric differential equations, called the basic Schwarz list. We refer to the mentioned 
for more details. 

We finish by a remark but first we recall well known facts on modular curves : the modular 
function J : fj t— > C gives a quotient map with respect to the projective group PSL(2,M). It 
ramifies above the points 0, 1 with ramification indices 3 and 2 respectively. The group T(m) — 
{M £ SL(2,Z), M = Id mod m} is a normal subgroup of SL(2,Z) and has no elliptic elements 
for to > 2. If Y(m) denotes the quotient of Sj by T(to), the cover f) H ► Y(m) is unramified and the 
modular function J factors over Y(m), J : Y(m) i— > C. The modular curve X(m) is the completion 
of Y(m) by adding to Y(m) cusps so that J : X(m) H> P 1 . This map ramifies above oo with order 
to. Hence the ramifications indices above 0, 1, oo are 3, 2, to. The covering group is PSL(2, Z/mZ). 
When to = 3, 4, 5 we recover the tetrahedral, octahedral and icosahedral coverings. 

In the previous theorem, the vanishing of the fourth transvectant appeared as one of the con- 
ditions for a solution of a Riccati equation to be algebraic. In the analysis of this vanishing, the 
next theorem is of great importance and it is due to Brioschi and Wedekind. 

Theorem 28 (0,01). Let K be an algebraically closed field of characteristic zero. The fourth 
transvectant (/, J)*- 4 - 1 of a non-trivial binary form f of order k > 4 is identically zero if and only if 
f is GL(2, K)-equivalent to one of the following forms 

(1) x\ orx\~ 1 X2 (degenerate case) 

(2) X2{x\ + x%) (tetrahedral case) 

(3) XiX2 {x\ + xf) (octahedral case) 

(4) X1X2 (x\° — llxfx! — ^2°) (icosahedral case). 



The main point here is that the vanishing of the fourth transvectant of a non-degenerate binary 
form / forces it to be one of the Klein forms. 
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6. Note 

There are various way one could possibly extend the study realized in this paper. First we may 
consider the case of homogeneous linear third order or n-th order 

S/ (n) + a lV ( n -V + . . . + a n y = 0, 

y> 

n > 2 over (C(x),S), 5(x) = 1 and C algebraically closed. Using the substitution u — — one reduces 

the homogeneous linear differential equation to a non-linear differential equation of order n — 1 . 
Do there exist anharmonics in this case? Moreover does the Darboux formalism in the last section 
of this paper generalize to the case of Riccati equations on the field of functions of an algebraic 
curve over C (the case we studied is essentially the case of elliptic curves). Finally one possible 
further generalization of the Darboux formalism may occur through the study of the general Abel 
differential equations 

du 



with q in some function field. 



T = u ' 3 + q 
dz 
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